function mpc = pglib_opf_case39_epri
mpc.version = '2';
mpc.baseMVA = 100;

%% bus data
%	bus_i	type	Pd	Qd	Gs	Bs	area	Vm	Va	baseKV	zone	Vmax	Vmin
mpc.bus = [
	1	 1	 97.6	 44.2	 0.0	 0.0	 2	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
	2	 1	 0.0	 0.0	 0.0	 0.0	 2	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
	3	 1	 322.0	 2.4	 0.0	 0.0	 2	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
	4	 1	 1000.0	 368.0	 0.0	 0.0	 1	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
	5	 1	 0.0	 0.0	 0.0	 0.0	 1	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
	6	 1	 0.0	 0.0	 0.0	 0.0	 1	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
	7	 1	 233.8	 84.0	 0.0	 0.0	 1	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
	8	 1	 1044.0	 353.2	 0.0	 0.0	 1	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
	9	 1	 6.5	 -66.6	 0.0	 0.0	 1	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
	10	 1	 0.0	 0.0	 0.0	 0.0	 1	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
	11	 1	 0.0	 0.0	 0.0	 0.0	 1	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
	12	 1	 8.53	 88.0	 0.0	 0.0	 1	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
	13	 1	 0.0	 0.0	 0.0	 0.0	 1	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
	14	 1	 0.0	 0.0	 0.0	 0.0	 1	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
	15	 1	 320.0	 153.0	 0.0	 0.0	 3	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
	16	 1	 658.0	 64.6	 0.0	 0.0	 3	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
	17	 1	 0.0	 0.0	 0.0	 0.0	 2	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
	18	 1	 158.0	 30.0	 0.0	 0.0	 2	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
	19	 1	 0.0	 0.0	 0.0	 0.0	 3	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
	20	 1	 680.0	 103.0	 0.0	 0.0	 3	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
	21	 1	 274.0	 115.0	 0.0	 0.0	 3	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
	22	 1	 0.0	 0.0	 0.0	 0.0	 3	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
	23	 1	 247.5	 84.6	 0.0	 0.0	 3	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
	24	 1	 1360.6	 206	 0.0	 0.0	 3	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
	25	 1	 224.0	 47.2	 0.0	 0.0	 2	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
	26	 1	 139.0	 17.0	 0.0	 0.0	 2	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
	27	 1	 281.0	 75.5	 0.0	 0.0	 2	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
	28	 1	 206.0	 27.6	 0.0	 0.0	 3	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
	29	 1	 283.5	 26.9	 0.0	 0.0	 3	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
	30	 2	 0.0	 0.0	 0.0	 0.0	 2	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
	31	 3	 9.2	 4.6	 0.0	 0.0	 1	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
	32	 2	 0.0	 0.0	 0.0	 0.0	 1	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
	33	 2	 0.0	 0.0	 0.0	 0.0	 3	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
	34	 2	 0.0	 0.0	 0.0	 0.0	 3	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
	35	 2	 0.0	 0.0	 0.0	 0.0	 3	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
	36	 2	 0.0	 0.0	 0.0	 0.0	 3	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
	37	 2	 0.0	 0.0	 0.0	 0.0	 2	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
	38	 2	 0.0	 0.0	 0.0	 0.0	 3	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
	39	 2	 1104.0	 250.0	 0.0	 0.0	 1	    1.00000	    0.00000	 345.0	 1	    1.06000	    0.94000;
];

%% generator data
%	bus	Pg	Qg	Qmax	Qmin	Vg	mBase	status	Pmax	Pmin	Pc1	Pc2	Qc1min	Qc1max	Qc2min	Qc2max	ramp_agc	ramp_10	ramp_30	ramp_q	apf
mpc.gen = [
	30	 520.0	 270.0	 800.0	  -800.0	 1.0499	 100.0	 1	 2080.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0; % NUC
	31	 323.0	 100.0	 300.0	 -100.0	 0.982	 100.0	 1	 646.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0; % COW
	32	 362.5	 225.0	 600.0	  -600.0	 0.9841	 100.0	 1	 1450.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0; % COW
	33	 326.0	 125.0	 250.0	  0.0	   0.9972	 100.0	 1	 652.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0; % COW
	34	 254.0	 83.5		 167.0	  0.0	   1.0123	 100.0	 1	 508.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0; % COW
	35	 343.5	 100.0	 600.0	 -600.0	 1.0494	 100.0	 1	 1374.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0; % COW
	36	 290.0	 120.0	 240.0	  0.0	   1.0636	 100.0	 1	 580.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0; % COW
	37	 282.0	 125.0	 250.0	  0.0	   1.0275	 100.0	 1	 564.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0; % COW
	38	 432.5	 75.0		 600.0	 -600.0	 1.0265	 100.0	 1	 1730.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0; % COW
	39	 550.0	 100.0	 300.0	 -100.0	 1.03	   100.0	 1	 1100.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0	 0.0; % COW
];



%% branch data
%	fbus	tbus	r	x	b	rateA	rateB	rateC	ratio	angle	status	angmin	angmax
mpc.branch = [
	1	 2	 0.0035	 0.0411	 0.6987	 600.0	 600.0	 600.0	 0.0	 0.0	 1	 -30.0	 30.0;
	1	 39	 0.001	 0.025	 0.75	 1000.0	 1000.0	 1000.0	 0.0	 0.0	 1	 -30.0	 30.0;
	2	 3	 0.0013	 0.0151	 0.2572	 500.0	 500.0	 500.0	 0.0	 0.0	 1	 -30.0	 30.0;
	2	 25	 0.007	 0.0086	 0.146	 500.0	 500.0	 500.0	 0.0	 0.0	 1	 -30.0	 30.0;
	2	 30	 0.0	 0.0181	 0.0	 900.0	 900.0	 2500.0	 1.025	 0.0	 1	 -30.0	 30.0;
	3	 4	 0.0013	 0.0213	 0.2214	 500.0	 500.0	 500.0	 0.0	 0.0	 1	 -30.0	 30.0;
	3	 18	 0.0011	 0.0133	 0.2138	 500.0	 500.0	 500.0	 0.0	 0.0	 1	 -30.0	 30.0;
	4	 5	 0.0008	 0.0128	 0.1342	 600.0	 600.0	 600.0	 0.0	 0.0	 1	 -30.0	 30.0;
	4	 14	 0.0008	 0.0129	 0.1382	 500.0	 500.0	 500.0	 0.0	 0.0	 1	 -30.0	 30.0;
	5	 6	 0.0002	 0.0026	 0.0434	 1200.0	 1200.0	 1200.0	 0.0	 0.0	 1	 -30.0	 30.0;
	5	 8	 0.0008	 0.0112	 0.1476	 900.0	 900.0	 900.0	 0.0	 0.0	 1	 -30.0	 30.0;
	6	 7	 0.0006	 0.0092	 0.113	 900.0	 900.0	 900.0	 0.0	 0.0	 1	 -30.0	 30.0;
	6	 11	 0.0007	 0.0082	 0.1389	 480.0	 480.0	 480.0	 0.0	 0.0	 1	 -30.0	 30.0;
	6	 31	 0.0	 0.025	 0.0	 1800.0	 1800.0	 1800.0	 1.07	 0.0	 1	 -30.0	 30.0;
	7	 8	 0.0004	 0.0046	 0.078	 900.0	 900.0	 900.0	 0.0	 0.0	 1	 -30.0	 30.0;
	8	 9	 0.0023	 0.0363	 0.3804	 900.0	 900.0	 900.0	 0.0	 0.0	 1	 -30.0	 30.0;
	9	 39	 0.001	 0.025	 1.2	 900.0	 900.0	 900.0	 0.0	 0.0	 1	 -30.0	 30.0;
	10	 11	 0.0004	 0.0043	 0.0729	 600.0	 600.0	 600.0	 0.0	 0.0	 1	 -30.0	 30.0;
	10	 13	 0.0004	 0.0043	 0.0729	 600.0	 600.0	 600.0	 0.0	 0.0	 1	 -30.0	 30.0;
	10	 32	 0.0	 0.02	 0.0	 900.0	 900.0	 2500.0	 1.07	 0.0	 1	 -30.0	 30.0;
	12	 11	 0.0016	 0.0435	 0.0	 500.0	 500.0	 500.0	 1.006	 0.0	 1	 -30.0	 30.0;
	12	 13	 0.0016	 0.0435	 0.0	 500.0	 500.0	 500.0	 1.006	 0.0	 1	 -30.0	 30.0;
	13	 14	 0.0009	 0.0101	 0.1723	 600.0	 600.0	 600.0	 0.0	 0.0	 1	 -30.0	 30.0;
	14	 15	 0.0018	 0.0217	 0.366	 600.0	 600.0	 600.0	 0.0	 0.0	 1	 -30.0	 30.0;
	15	 16	 0.0009	 0.0094	 0.171	 600.0	 600.0	 600.0	 0.0	 0.0	 1	 -30.0	 30.0;
	16	 17	 0.0007	 0.0089	 0.1342	 600.0	 600.0	 600.0	 0.0	 0.0	 1	 -30.0	 30.0;
	16	 19	 0.0016	 0.0195	 0.304	 600.0	 600.0	 2500.0	 0.0	 0.0	 1	 -30.0	 30.0;
	16	 21	 0.0008	 0.0135	 0.2548	 600.0	 600.0	 600.0	 0.0	 0.0	 1	 -30.0	 30.0;
	16	 24	 0.0003	 0.0059	 0.068	 600.0	 600.0	 600.0	 0.0	 0.0	 1	 -30.0	 30.0;
	17	 18	 0.0007	 0.0082	 0.1319	 600.0	 600.0	 600.0	 0.0	 0.0	 1	 -30.0	 30.0;
	17	 27	 0.0013	 0.0173	 0.3216	 600.0	 600.0	 600.0	 0.0	 0.0	 1	 -30.0	 30.0;
	19	 20	 0.0007	 0.0138	 0.0	 900.0	 900.0	 2500.0	 1.06	 0.0	 1	 -30.0	 30.0;
	19	 33	 0.0007	 0.0142	 0.0	 900.0	 900.0	 2500.0	 1.07	 0.0	 1	 -30.0	 30.0;
	20	 34	 0.0009	 0.018	 0.0	 900.0	 900.0	 2500.0	 1.009	 0.0	 1	 -30.0	 30.0;
	21	 22	 0.0008	 0.014	 0.2565	 900.0	 900.0	 900.0	 0.0	 0.0	 1	 -30.0	 30.0;
	22	 23	 0.0006	 0.0096	 0.1846	 600.0	 600.0	 600.0	 0.0	 0.0	 1	 -30.0	 30.0;
	22	 35	 0.0	 0.0143	 0.0	 900.0	 900.0	 2500.0	 1.025	 0.0	 1	 -30.0	 30.0;
	23	 24	 0.0022	 0.035	 0.361	 600.0	 600.0	 600.0	 0.0	 0.0	 1	 -30.0	 30.0;
	23	 36	 0.0005	 0.0272	 0.0	 900.0	 900.0	 2500.0	 0.0	 0.0	 1	 -30.0	 30.0;
	25	 26	 0.0032	 0.0323	 0.531	 600.0	 600.0	 600.0	 0.0	 0.0	 1	 -30.0	 30.0;
	25	 37	 0.0006	 0.0232	 0.0	 900.0	 900.0	 2500.0	 1.025	 0.0	 1	 -30.0	 30.0;
	26	 27	 0.0014	 0.0147	 0.2396	 600.0	 600.0	 600.0	 0.0	 0.0	 1	 -30.0	 30.0;
	26	 28	 0.0043	 0.0474	 0.7802	 600.0	 600.0	 600.0	 0.0	 0.0	 1	 -30.0	 30.0;
	26	 29	 0.0057	 0.0625	 1.029	 600.0	 600.0	 600.0	 0.0	 0.0	 1	 -30.0	 30.0;
	28	 29	 0.0014	 0.0151	 0.249	 600.0	 600.0	 600.0	 0.0	 0.0	 1	 -30.0	 30.0;
	29	 38	 0.0008	 0.0156	 0.0	 1200.0	 1200.0	 2500.0	 1.025	 0.0	 1	 -30.0	 30.0;
];

%% candidate dc bus data
%column_names%   busdc_i grid    Pdc     Vdc     basekVdc    Vdcmax  Vdcmin  Cdc
mpc.busdc_ne = [
		1              1       0       1       345         1.1     0.9     0;
		2              1       0       1       345         1.1     0.9     0;
		3              1       0       1       345         1.1     0.9     0;
		4              1       0       1       345         1.1     0.9     0;
		5              1       0       1       345         1.1     0.9     0;
		6              1       0       1       345         1.1     0.9     0;
		7              1       0       1       345         1.1     0.9     0;
		8              1       0       1       345         1.1     0.9     0;
];

%% candidate branches
%column_names%   fbusdc  tbusdc  r      l        c   rateA   rateB   rateC status cost
mpc.branchdc_ne = [
1	5      0.01	 0.00	 0.00  400.0	 0.0	 0.0	 1.0	1.1;
1	6      0.01	 0.00	 0.00  400.0	 0.0	 0.0	 1.0	1.2;
1	7      0.01	 0.00	 0.00  400.0	 0.0	 0.0	 1.0	1.3;
1	8      0.01	 0.00	 0.00  400.0	 0.0	 0.0	 1.0	1.4;
2	5      0.01	 0.00	 0.00  400.0	 0.0	 0.0	 1.0	1.5;
2	6      0.01	 0.00	 0.00  400.0	 0.0	 0.0	 1.0	1.6;
2	7      0.01	 0.00	 0.00  400.0	 0.0	 0.0	 1.0	1.7;
2	8      0.01	 0.00	 0.00  400.0	 0.0	 0.0	 1.0	1.8;
3	5      0.01	 0.00	 0.00  400.0	 0.0	 0.0	 1.0	1.9;
3	6      0.01	 0.00	 0.00  400.0	 0.0	 0.0	 1.0	2;
3	7      0.01	 0.00	 0.00  400.0	 0.0	 0.0	 1.0	2.1;
3	8      0.01	 0.00	 0.00  400.0	 0.0	 0.0	 1.0	2.2;
4  5      0.01	 0.00	 0.00  400.0	 0.0	 0.0	 1.0	2.3;
4  6      0.01	 0.00	 0.00  400.0	 0.0	 0.0	 1.0	2.4;
4	7      0.01	 0.00	 0.00  400.0	 0.0	 0.0	 1.0	2.5;
4  8      0.01	 0.00	 0.00  400.0	 0.0	 0.0	 1.0	2.6;

1	5       0.01	 0.00	 0.00  800.0	 0.0	 0.0	 1.0	3.1;
1	6       0.01	 0.00	 0.00  800.0	 0.0	 0.0	 1.0	3.2;
1	7       0.01	 0.00	 0.00  800.0	 0.0	 0.0	 1.0	3.3;
1	8       0.01	 0.00	 0.00  800.0	 0.0	 0.0	 1.0	3.4;
2	5       0.01	 0.00	 0.00  800.0	 0.0	 0.0	 1.0	3.5;
2	6       0.01	 0.00	 0.00  800.0	 0.0	 0.0	 1.0	3.6;
2	7       0.01	 0.00	 0.00  800.0	 0.0	 0.0	 1.0	3.7;
2	8       0.01	 0.00	 0.00  800.0	 0.0	 0.0	 1.0	3.8;
3	5       0.01	 0.00	 0.00  800.0	 0.0	 0.0	 1.0	3.9;
3	6       0.01	 0.00	 0.00  800.0	 0.0	 0.0	 1.0	4;
3	7       0.01	 0.00	 0.00  800.0	 0.0	 0.0	 1.0	4.1;
3	8       0.01	 0.00	 0.00  800.0	 0.0	 0.0	 1.0	4.2;
4  5        0.01	 0.00	 0.00  800.0	 0.0	 0.0	 1.0	4.3;
4  6        0.01	 0.00	 0.00  800.0	 0.0	 0.0	 1.0	4.4;
4	7       0.01	 0.00	 0.00  800.0	 0.0	 0.0	 1.0	4.5;
4  8        0.01	 0.00	 0.00  800.0	 0.0	 0.0	 1.0	4.6;

1	5       0.01	 0.00	 0.00  1200.0	 0.0	 0.0	 1.0	5.1;
1	6       0.01	 0.00	 0.00  1200.0	 0.0	 0.0	 1.0	5.2;
1	7       0.01	 0.00	 0.00  1200.0	 0.0	 0.0	 1.0	5.3;
1	8       0.01	 0.00	 0.00  1200.0	 0.0	 0.0	 1.0	5.4;
2	5       0.01	 0.00	 0.00  1200.0	 0.0	 0.0	 1.0	5.5;
2	6       0.01	 0.00	 0.00  1200.0	 0.0	 0.0	 1.0	5.6;
2	7       0.01	 0.00	 0.00  1200.0	 0.0	 0.0	 1.0	5.7;
2	8       0.01	 0.00	 0.00  1200.0	 0.0	 0.0	 1.0	5.8;
3	5       0.01	 0.00	 0.00  1200.0	 0.0	 0.0	 1.0	5.9;
3	6       0.01	 0.00	 0.00  1200.0	 0.0	 0.0	 1.0	6;
3	7       0.01	 0.00	 0.00  1200.0	 0.0	 0.0	 1.0	6.1;
3	8       0.01	 0.00	 0.00  1200.0	 0.0	 0.0	 1.0	6.2;
4  5        0.01	 0.00	 0.00  1200.0	 0.0	 0.0	 1.0	6.3;
4  6        0.01	 0.00	 0.00  1200.0	 0.0	 0.0	 1.0	6.4;
4 7       0.01	 0.00	 0.00  1200.0	 0.0	 0.0	 1.0	6.5;
4  8        0.01	 0.00	 0.00  1200.0	 0.0	 0.0	 1.0	6.6;
];

%% candidate converters
%column_names%   busdc_i busac_i type_dc type_ac P_g   Q_g  islcc  Vtar rtf xtf  transformer tm   bf filter    rc      xc  reactor   basekVac Vmmax   Vmmin   Imax    status   LossA LossB  LossCrec LossCinv  droop Pdcset    Vdcset  dVdcset Pacmax Pacmin Qacmax Qacmin cost
mpc.convdc_ne = [
	1       4   2       1       -360    -1.66   		0 1.0        0.01  0.01 1 1 0.01 1 0.01   0.01 1  345         1.1     0.9     15     1      1.1033 0.887  2.885    2.885        0.0050    -52.7   1.0079   0  2500 -2500 2500 -2500		3.1;
	2       8   2       1       -360    -1.66   		0 1.0        0.01  0.01 1 1 0.01 1 0.01   0.01 1  345         1.1     0.9     15     1      1.1033 0.887  2.885    2.885         0.0050    -52.7   1.0079   0  2500 -2500 2500 -2500	3.2;
	3      16   2       1       -360    -1.66   		0 1.0        0.01  0.01 1 1 0.01 1 0.01   0.01 1  345         1.1     0.9     15     1      1.1033 0.887  2.885    2.885        0.0050    -52.7   1.0079   0  2500 -2500 2500 -2500		3.3;
	4      24   2       1       -360    -1.66   		0 1.0        0.01  0.01 1 1 0.01 1 0.01   0.01 1  345         1.1     0.9     15     1      1.1033 0.887  2.885    2.885        0.0050    -52.7   1.0079   0  2500 -2500 2500 -2500		3.4;
	5      30   2       1       -360    -1.66   		0 1.0        0.01  0.01 1 1 0.01 1 0.01   0.01 1  345         1.1     0.9     15     1      1.1033 0.887  2.885    2.885        0.0050    -52.7   1.0079   0  2500 -2500 2500 -2500		3.5;
	6      32   2       1       -360    -1.66   		0 1.0        0.01  0.01 1 1 0.01 1 0.01   0.01 1  345         1.1     0.9     15     1      1.1033 0.887  2.885    2.885        0.0050    -52.7   1.0079   0  2500 -2500 2500 -2500		3.6;
	7    	 35   2       1       -360    -1.66   		0 1.0        0.01  0.01 1 1 0.01 1 0.01   0.01 1  345         1.1     0.9     15     1      1.1033 0.887  2.885    2.885         0.0050    -52.7   1.0079   0  2500 -2500 2500 -2500	3.7;
	8      38   2       1       -360    -1.66   		0 1.0        0.01  0.01 1 1 0.01 1 0.01   0.01 1  345         1.1     0.9     15     1      1.1033 0.887  2.885    2.885        0.0050    -52.7   1.0079   0  2500 -2500 2500 -2500		3.8;
	];

%% generator cost data
%	2	startup	shutdown	n	c(n-1)	...	c0
mpc.gencost = [
	2	 0.0	 0.0	 3	   0.000000	  0.001   0.000000; % NUC
	2	 0.0	 0.0	 3	   0.000000	  0.001   0.000000; % COW
	2	 0.0	 0.0	 3	   0.000000	  0.001   0.000000; % COW
	2	 0.0	 0.0	 3	   0.000000	  0.001   0.000000; % COW
	2	 0.0	 0.0	 3	   0.000000	  0.001   0.000000; % COW
	2	 0.0	 0.0	 3	   0.000000	  0.001   0.000000; % COW
	2	 0.0	 0.0	 3	   0.000000	  0.001   0.000000; % COW
	2	 0.0	 0.0	 3	   0.000000	  0.001   0.000000; % COW
	2	 0.0	 0.0	 3	   0.000000	  0.001   0.000000; % COW
	2	 0.0	 0.0	 3	   0.000000	  0.001   0.000000; % COW
];
